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abstract. De riaan ana fereira 1 2006 ) provided models for spatial extremes in the case of stationarity, 



which depend on just one parameter > measuring tail dependence, and they proposed different estimators 



for this parameter. This framework was supplemented in IFalkl 1 2011 ) by establishing local asymptotic 
normality (LAN) of a corresponding point process of exceedances above a high multivariate threshold, yielding 
in particular asymptotic efficient estimators. 

The estimators investigated in these papers are based on a finite set of points t\ , . . . , tj, at which observa- 
tions are taken. We generalize this approach in the context of functional extreme value theory (EVT). This 
more general framework allows estimation over some spatial parameter space, i.e., the finite set of points 
ti, . . . ,td is replaced by t G [a,b]. In particular, we derive efficient estimators of /3 based on those processes 
in a sample of iid processes in C[0, 1] which exceed a given threshold function. 



1. Introduction 



Supp ose that the stochasti c process V = (Vf) te [ ,i] € C[0, 1] is a standard generalized Pareto process 
(GPP) (jBuishand et al.l (|2008l )). i.e., there exists xq > such that 



P(V</) = l + log(G(/)), feE-[0,l], H/IL^zo, 

where E~[0, 1] is the set of those bounded functions on [0, 1] that attain only nonpositive values and which 
have a finite set of discontinuities. By G we denote the functional distribution function (df) of a standard 
extreme value process (EVP) rj = (?7t)te[o.i] S C[0, 1], i.e., 

G(f) = P( V <f), f€E-[0,l], 
P(i] t < x) = cxp(x), x < 0, t G [0,1], and rj is max-stable: 

?(*)<£) =P(v<f), /e£-[o,i],neN. 

All operations on functions such as <, multiplication with a constant etc. are meant componentwise. For 
random functions, i.e., stochastic processes such as V, r\ we use bold letters, to distinguish these from 
nonrandom functions such as f. 

De Haan and Pereira (|2006l ) provided models for spatial extremes in the case of stationarity, which depend 
on just one parameter j3 > measuring tail depe ndence, and they proposed different estimators for this 
parameter. This framework was supplemented in IFalkl (|201ll ) by establishing local asymptotic normality 
(LAN) of a corresponding point process of exceedances above a high multivariate threshold. 
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Precisely, it is assumed that for any x\ 



,x d < 0, d G N, 



(1) 



PiVtj < Xj, 1 < j < d) = exp I - 



max \ x, 

3<d 



\ipp{s - tj)ds 



where ipp (s) = /3ip(/3s) with a scale parameter j3 > 0, and tp is a continuous probability density on 
?/>(s) = ip{—s) > a nd ^(a), s > 0, decreasing. 



with 



In the papers by Ide Haan and Pereiral (|2006[ ) and iFalkl (|201lf ) the density ip is known and the parameter 
f3 is estimated. The estimators investigated in these papers are based on a finite set of points t\ < ■ ■ ■ < td] 
estimation over some interval t G [a, b] seems to be an open problem. This is the content of the present 
paper, which is organized as follows. In Section [5] wc compile some auxiliary results and tools, in particular 
from functional extreme value theory (EVT). In Section [3] we introduce our estimator of /3 and establish its 
asymptotic normality under the condition that the underlying observations V^', . . . , V( n > are independent 
copies of a standard GPP V. Local asymptotic normality (LAN) of a corresponding point process of ex- 
ceedances above a high constant threshold function is established in Section SJ This is achieved under the 
condition that the underlying observations are in a (5-neighborhood of a standard GPP. As an application we 
obtain from LAN-theory th at our estimator of P i s asym ptotically efficient in this setup. For an account of 
functional EVT we refer to 
GPP we refer to 



Aulbach et al 



poiih . 



de Haan and Ferreiral (|2006f ); for a supplement including in particular basics of 



2. Auxiliary Results and Tools 

In this section we compile several auxiliary results and tools. We start with the functional df of a standard 
EVP r] G C[0, 1], whose finite dimensional marginal distributions (fidis) are given by equation ([T]). 



Lemma 2.1. We have for any f G E [0, 1] 



P{ri <f) = cxp (- J 



sup (\f(t)\^(s- pt)) ds 
oote[o,i] , 



Proof. The assertion follows from the fact that a probability measure is continuous from above together with 
the dominated convergence theorem; note that J_ oa sup te [ y ip(s — fit) ds < oo. Let Q = {91,92, ■ • • } be 
a denumerable and dense subset of [0,1], which contains also the set of discontinuities of /. Recall that 
T] G C[0, 1]. From representation (fTJ) we obtain 



p(v < /) = p ( n ^ /(»)> 1 < * < ^ ) 

VriGN / 

= lim P(r) qi < /(%), 1 < i < n) 

n— too 



lim exp 



exp 



exp 



max (\f(qi)\ tpp(s - qi)) ds 

Ki<n 



lim max (|/(<? 4 )| ippi-s - q t )) ds 

n— foo \ Ki<n 



sup (\f(t)\^(s-f3t)) ds 

-ootG[0,l] , 



□ 



The preceding result provides the functional df P(V </) = ! + log(G(/)) of the GPP V in its upper tail. 
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Corollary 2.2. There exists xq > such that for the GPP V corresponding to the EVP r\ and for any 
feE~[0,l] with ll/H^ < x 

(i) 

/oo 
sup (\f(t)\i>(s-0t)) ds, 
-oo te[o,il 



(ii) 



/OO 
inf (\f(t)\il>(s-0t)) ds. 
-oo*e[o,i] 



Proof. While part (i) is an immediate consequen ce of of Lem ma 12.11 part (ii) follows from the inclusion- 
exclusion formula as in the proof of Lemma 3.1 in iFalk ( 2011 '). □ 



Note that 



/CO 
sup (\f(t)\i/>(s-pf)) ds 7 
-oo te[o,il 



te[o,i] 

defines a norm on the set E[0, 1], called D-norm. By E[0, 1] we denote the set of those functions on [0, 1], 
which are bounded and have a finite number of discontinuities. The representation of a multivariate extreme 
value distributi on (EVP) or of a multivariate generalized Pareto distribution (GPD) i n terms of a .D-norm i s 
well-known, see 



Falk et al 



(|2010t ) . This concept was extended to functional spaces in lAulbach et al.l (|201lf ). 
The fidis of the stochastic processes r/orV are obtained by considering the function f(t) = £^ =1 Xilr t .y(t) G 
^"[O, 1], Xi < 0, ti e [0, 1], d € N. This norm satisfies, for example, the general inequality 



ll/IL< II/IId 

where 1 denotes the constant function one and 



< 



\D > 



feE[0,l], 



x := sup te [ ^ \ f{t)\ is the usual sup-norm. This inequality 
implies in particular that each D-norm is equivalent with the sup-norm which, in turn, implies that the L p - 

norm ||/|| = (J \.f(t)\ p dt\ , with p g [l,oo), is not a D-norm. 

The following auxiliary result is a crucial tool for the derivation of estimators of f3. 



Lemma 2.3. We have 

where ^/(x) = J_ ip( s ) ds. 
Proof. We have 



inf ip(s — fit) ds 
te[o,i] 



inf ibis — fit) ds 
oo*e[o,i] 



2 1 - # - = 2* -- 



inf ib(\s - Bt\) ds 
oo*e[o,i] 

so 

min (^(|s|), ip{\s — (3\)) ds 



ip{s) ds ■ 

0/2 

2*(-0/2). 



0/2 



ip(s-0) ds 



□ 
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3. Estimation of p 

A natural estimator of P/2), based on independent copies V«, . . . , VW of V, is by Corollary [22 
and Lemma 12.31 given by 



Note the twofold meaning of c: In the denominator 2 |c| n this is just the absolute value of the constant c < 0, 
whereas in the term l( c .o](V r< - I - ) ) we mean the constant function c, and we have l( c ,o](V''') = 1 if and only if 
each component satisfies v} 1 ^ > c, t G [0, 1]. There should be no risk of confusion. 
The law of large number implies 

$c,n -> n ->oo * a.S. 

and, thus, 

/3 c>n := -2*" 1 f$ c J ^„_>oc, P a.s., 

where F^ 1 (q) := inf {fel: F(t) > q}, q G (0, 1), denotes the generalized inverse of a df F. 

The Moivre-Laplace theorem implies as ympt o tic no rmality of ^ c ^ n and P c .m i.e., the next result is a 



functional counterpart of Proposition 3.3 in iFalkl (|201l[ ). 
Proposition 3.1. For c < close enough to we have 



r D 



N 0, 



ttf-f l-2|c|* -f 



and 



l 1/2 (Pc,n ~ P) N I 0, 



2 Id 



-f j (l-2|c|*(-f 



|d> 2 (-1) 



We now consider a stochastic process X G (^"[O,!] := {/ G C[0, 1] : / < 0}, whose upper tail is in 
a ^-neighborhood of that of a GPP V G C~[0,1] with Z?-norm \\f\\ D = sup te [ 0il ](|/(t)| ip(s — pt))ds. 
Precisely, we require that 

(C) P(X > cf) = P(V > cf) (1 + c s K(f) + r(c, /)) 



for c G (0, 1) and / G E [0, 1] with ||/|| < £o for some e > 0, where K : E [0, 1] — > M is a function and 
the remainder r(c, /) is of order o (c 5 ) as c — > 0. The next result is an immediate consequence of Corollary 
2~2l 



Lemma 3.2. Suppose that the stochastic process X G C [0, 1] satisfies condition ((C|). Then we obtain for 
cg (0,1) andfEE-[0,l] with WfW^ < e 

P(X > cf) =c( f°° inf (\f(t)\ Ms - pt)) ds) (1 + c s K(f) + r(c, /)) . 
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In w hat follows we show how a process X satisfying condition ([C]) can be generated. From I Aulbach et al 
(|201l[ ) we conclude that there is a stochastic process Z = (Z t ) tg [ .i] on [0, 1] with continuous sample paths 
and < Z t < m, E(Z t ) = 1, t G [0, 1], for some constant m > 1, such that 



\\f\\ D = / sup (\f(t)\^(s-f3t))ds 
J-oo te[o,i] 

= E[ sup (\f(t)\Z t ) ) , f€E[0,l). 
\te[o,i] J 

The stochastic process Z is called generator of the D-norm. Conversely, each process Z with the above 
properties generates a D-norm via \\,f\\ D := E ^sup te [ 01 ] (|/(i)| Z t )j , / G E[0, 1]. For every D-norm 
there exists a standard EVP 77 e C[0, 1] with functional df P(r] < /) = exp(- \\f\\ D ), f G ^"[0, 1]. While a 
generat or Z is in general not uniquely determined, the generator constant E ^sup tg j = \\1\\ D is. We 

refer to 
Put 



Aulbach et al 



(|201lh for details. 



(2) V := (y t ) tg [ 0l i] := fmaxf~M 

V V z * //te[o,i] 

where [/ and ^ arc independent, U is a uniformly on (0, 1) distributed rv and M < is an arbitrary constant. 
We incorporate the constant M to ensure that V t > —00 for each t G [0, 1], as Z t may attain the value zero. 
The continuous process V is a GPP, as we have for / G E~[0, 1] with WfW^ < min(|M|, 1/m) 

^P(V < /) = P(f/ > |/(f)| Z t ,Q<t<i) 

= P(U> sup (|/(t)|Z t )) 
te[o,i] 

= l-^[ sup (|/(t)|^t) J 

\*6[o,i] y 

= 1-11/112,- 

We have, moreover, 

P(V>/) = P(£7<|/(t)|Z t , < * < 1) 



E inf (|/(t)|Z t ) 
Ve[o,i] 

inf (|/(t)|^( S -^))ds, 

-ootGtO.l] 

where the final equality is a consequence of Corollarv l2.2i part (ii). 

Replace now the rv U in ([2]) by a rv F > 0, which is also independent of Z and whose df H is continuous 
and satisfies 

(3) H{u) =u + Au 1+S + o(u 1+s ) as u I 

with some constant The standard exponential distribution, for instance, satisfies this condition with 

5 = 1 and A = —1/2. The process 

Y 

7j t //te[o,i] 



(4) X := (X t ) te[0 ,i] := ( max ( -— ,M 
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then satisfies condition ((C|) with 

K( f) = A E(M tmi] (\f(t)\Z t )^) 
E(M te[0A] (\f(t)\Z t )) ' 

which has to be interpreted as zero if the denominator vanishes. 

The following theorem is the main result of this section. We will see in Section 2] using LAN theory that 
it implies that & Cn n is an asymptotically efficient estimator sequence in an appropriate model. 



Theorem 3.3. Suppose that the stochastic process XgC [0, 1] satisfies condition fCj. // the sequence of 
hoh 

(i) 



thresholds c n < 0, n £ N, satisfies c n — > ; n \c n \ — > oo, n c„| 1+2<5 — > const > as n — > oo, then we obtain 



(n M) 1/2 (*c,» - * (-§ )) N (constVV ^ (-§)) , 

(ii) 

, , iM/2/3 o\ „( Sconst 1 /^ 2*(-|)\ 

where fx := K(-T)V(-0/2). 
Proof. From Lemma 13.21 we obtain 



1 fP(X > 



cf V M 



c 



inf ip(s — fit) ds 



-+ cf0 K(-l) I t w£ 4(a - [3t) ds 



(5) = 2ff(-l)tt(-£ 



Write 



= (n\c n \) 1 - 

= -Vn+ K 

The Moivrc-Laplace theorem implies 



(4^g (i <-»i (Xj, - p(x>c " ,) ) 



and expansion ([5]) yields 



(n| 



cJ 1+25 



1/2 



1 /P(X>c„) 



6„ = t, ■ — i — i / inf ^(s - (3t) ds 



\cJ 5 V |Cn| J-oo*£[04] 



-^n-^-oo COnst 1 / 2 ^. 
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Equally, one concludes 



2(n| Cn |) 1 /2(v & -i(v & (-f))-vI,-i($ c 



2 

2(n|c„|) 1 / 2 (tt 

2const 1 /2 M 2*(-f 



— >D iV I — - 



by Slutsky's lemma, with £ between \& c „,n and ^(—j3/2). This completes the proof. □ 

The idea suggests itself to substitute the constant threshold c by a suitable threshold function / e E~[0, 1] 
and to consider, with c < 0, 

n 

11 z— 1 

^™ 2]^ P(F> |c|/) 
1 



inf (|/(t)|V(*-#)) ds 
2 J_oo*e[o,i] 



almost surely by the law of large numbers and Corollary [221 

The fact that with constant function / = — 1, the above integral equals by Lemma 



inf ib(s - Bt) ds = 2* , 

was the crucial observation for the derivation of an estimator of j3. Substituting the constant function / = — 1 
by an arbitrary function / <E E~ [0, 1] can, however, lead to surprising consequences, as the following example 
shows. 

Example 3.4. Take = exp(— \s\), sel, and f(t) := — exp(— t), t 6 [0,1]. Then we have for any 
B G [0, 1] 

inf (\f(t)\1>(8-pt)) ds = exp(-l), 

i.e., the functional T(f,f3,ip) is not capable to discriminate between different values of j3 G [0, 1]. For j3 > 1 
one obtains, however, 

1 + 



T(/,AV0 = exp 

The question, whether for each underlying density there exists an optimal threshold function / = /,/,, is 
an open problem. 

4. LAN OF EXCEEDANCES 

Let XW, 1 < i < n, be independent copies of a stochastic process X G ^"[0, 1], which satisfies condition 
(fCj) . Choose c < 0. In this section we establish local asymptotic normality (LAN) of the point process of 
exceedances 

n 

iV„, c (B):=^ £ , w n (fin [0,1)), Bel, 

* — ' su Pte[o,i] I A t / c j 
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where B denotes the cr-field of Borel sets in R. Note that for s 6 (0, 1] 

sup — < s •<=>■ X > sc, 
te[o,i] c 

i.e., the random point measure N n>c actually represents those processes among X^\ . . . , X^ which are 
exceedances above the constant function c. 

It is by Theorem 13.31 part (i), quite convenient to substitute the parameter /3 > in the family ippi') = 
f3?p(f3-) by the parameter 

^ := 2^ (-0 6(0,1). 

Fix $o 6 (0, 1). We require that the family of univariate df F#. c (s) := P®{X > sc), $ 6 (0, 1), s > 0, satisfies 
for s 6 (0, 1), Co < c < for some Co < 0, and # close to $o the expansion 

CD) : = K*?X >8C \ = 1 + ^ " **) + *, c), 

with 

r* o (M,c) = o(|i>-0o|) + O(|c| 7 ) 

uniformly for s 6 (0,1), Co < c < and $ close to i?o, where the constants L £ R and 7 > may 
depend on «?o- Note that condition (|D| implies in particular that P$ ^sup tg [ X ] (X t /c) = s^j = and, thus, 

^,c(s) = Pe ( su Pte[o,i] (^*/ c ) < s^) , s > 0, is actually a df on [0, 00). 

Condition fDJ is, for example, satisfied with L = l/i?o and r^ = if X is a GPP. We can also use the 
approach from definition (j4j) to generate a process 



X := ( max ( ~Y' M 



te[o,i] 



In addition to condition (j3|) we require that the continuous df H of the rv Y > satisfies the expansion 

H{u) = u + Au 1+s + r(u), 0<u<l, 

with some constant i e i, where the function r is diffcrentiable on (0, 1) with bounded derivative and 
r'(u) = o(u s ) as u 4. 0. Then condition fD]) is satisfied with L = l/#o- 

Denote by 11, ... , l^-( n ) those rv among sup tg [ ,i] /c, 1 < i < n, with sup tG r 0i i] /c < 1, in the order 
of their outcome. Then we have 

fe<r(n) 



From Theorem 1.4.1 in iReissI (|1993f ) we obtain that li, I2, . . . arc independent copies of a rv Y with df 

W<t) = ^^, 0<t<l, 
P&{X > c) 

under parameter 7? > 0, and they are independent of the total number r(n), which is binomial i?(n, P$(X > 
c))-distributed. 

Since the distribution C§ iC (Y) of V under -d is by condition ((D)) dominated by C^ 0iC (Y) for 1? in a neighbor- 
hood of $0 and cp < c < 0, the distribution C^(N niC ) of iV niC is dominated by (-/V„ jC ), see, e.g. Theorem 
3.1.2 in iReisd (|1993l ). Precisely, 7V„ ;C is a random element in the set M := {// = J2i<j< n £ Vj '■ Vi — 0> J — 



n, n = 0, 1, 2, . . . } of finite point measures on ([0, 00), B n [0, 00)), equipped with the smallest cr-field M. such 
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that for any B g B n JO, oo ) the projection ttb 



Section 1.1 in 
From 



{0, 1,2,...}, ttb(h) '■= n(B), is measurable; we refer to 



Reisa (|1993[ ) for technical details. 
Reisd (|1993l . Example 3.1.2) we conclude that C^(N n , c ) has the C$ (iV raiC )-density 
dCt(N n , c ) 



dC# (N n>c ) 
V(o,i)) 

n 



(/<) 



U,c(Vi) P Va( X > c ) 



P»jX > c) 
P# (X>c) 



M((0,1)) 



l-P d (X>c) 
l-P &0 {X>c) 



n-/i((0,l)) 



if (J, = Y^,i3i £ yi an( i /^((Oj 1)) < n - The loglikclihood ratio is, consequently, 



(t? I #o) 



log 



(JVn,c) 



fc<r(n) 



dCd {N n>c ) 

hJY k ) P^ (X>c) 



(7) 



r(n) lo, 



/* D ,c(n) p*(a:>c) 

P^(X > c 



+ (n — t(h)) log 



P* (X>c) 

We let in the sequel c = c„ depend on the sample size n with c„ 
as n . — > oo. Precisely, we put with arbitrary ^el 



>c) 

1-P, (X> C ) 
and, equally, d = i?„ with — > $o 



(n|c„|)V2' 



The following theorem is the main result of this section. It is analogous to Theorem 5.1 in lFaLa (|201lf ). whose 
proof carries over. 

Theorem 4.1. Suppose that ip(s) = ip(—s) and that ip(s), s > 0, is decreasing. Suppose, further, that 
n \c n \ — > n ^oo oo and that 

l+2min((5,7) 



(8) 

Then we obtain the expansion 



n c„ 



0. 



£n,c n (4 | t?o) = £LZ„ - 

-+°*0 N 



e^ 



eL 2 $o c2t2 



with 
(9) 



Z ™ — — KT7o >A) n iV (U, ^oj- 



(n|c„|)V2 

The above result reveals that the complete information about the underlying parameter that is contained 
in the exceedances Y\, . . . , Y T ( „\ is, actual ly, contained in their number r(n) as n increases. This is in complete 



accordance with the results in 



Falkl (|1998I ) , where this phenomenon was studied for general truncated empirical 



processes. The result here is, however, derived under more specialized conditions. 

Theorem 14. 1 1 together with the Hajek-LeCam convolution theorem provides the asymptotically minimum 
variance within the classes of regular estimators of $o- This class of estimators i?„ is defined by the property 
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that they are asymptotically unbiased under i} n = = $o + £,( n \ c n\) with i?o E (0, 1) for any 

precisely, 

(n |c„|) 1/2 (X - Qi> , 

where the limit distribution Q$ doe s not depend on £: see, e .g. Sections 8.4 and 8.5 in 

By LeCam's first lemma (see, e.g.. lLeCam and YansJ (|1990l . Chapter 3, Theorem 1)) we obtain that under 
0» = 0»(O 



Pfanzagl 



(|l994h . 



in,c„(4 I O )=£LZ n - 



£ 2 L 2 tfo 



with 
(10) 



. N(ZL# ,# ). 



An efficient estimator of $0 within the class of regular estimators has necessarily the minimum limiting 
variance 



2 

minimum 



1 



which is the inverse of the limiting variance of the central sequence LZ n under i?o (jPfanzagll (|1994 
8.4.1)). 

Consider the estimator 

r(n) 



Theorem 



n |c„| ' 



Then we have with t?„ = = tf + £( n |c„|)~ 1/2 



(n|c„|) 1/2 (? n -0 n ) =(«|c n |) 1 



/2 



n |c„| 



The estimator i) n is, consequently, not a regular estimator since we have by (|10p 

(n\c n \)^ 2 (d n -K) = Z n -£-> D , n N(t(M -l),d ), 

where the limiting distribution depends on £ unless L = l/f?o- Its asymptotic relative efficiency, defined as 
the ratio of the limiting variances under $0 is 

ARE(§o) = -2-^ 



minimum 



Recall that L = l/i?o if -X" follows a GPP or if X is in a neighborhood of a GPP as in (|6|) and, thus, i?„ is 
in this case regular and asymptotically efficient. 



Corollary 4.2. Suppose in addition to the conditions of Theorem ^. 1\ that X is a GPP or it is in a neigh- 
borhood of a GPP as in (JSJ). Then "& n = r(n)/(n \ c n \), n E N, is a regular estimator sequence with asymptotic 
minimum variance z?o within the class of regular estimators. 

A regular estimator sequence can in general be obtained as follows. Suppose that is a solution of the 
equation 

P K (X>c n ) = ^-. 
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Since t(h) is under z?o binomial B{n 1 P$ (X > c„)) distributed, $* is, actually, the maximum likelihood 
estimator of k = P$ {X > c n ) for the family {B(n, k) = B(n,P$(X > c„)) : i) £ (0,1)}. We suppose 
consistency of the sequence , n £ N. Then we obtain from condition (|D]l the expansion 

^ = p„. (x > cO 

n 

= / ( 1 + i(C _ ' !? o)+rtf (u,'!9*,c n ))/ l?0iCn (u)du 
Jo 

= (1 + L« - 0„) + O P , (\K ~M) + (|C„D) P» (X > Cn), 

which implies 

(n M) 1/2 - <?o) - + op, (1). 

As a consequence we obtain from © and (|10p with # n = i? n (£) 

(n |c„|) 1/2 (C - tf„) -> D , B A (o, ^2 ) . 

and, thus, $* , n G N, is a regular estimator sequence with asymptotic minimum variance. 
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